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In the solar convection zone, rotation couples with intensely turbulent convection to build global-scale flows of differential 
rotation and meridional circulation. Our sun must have rotated more rapidly in its past, as is suggested by observations 
of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the patterns of 
convection in such stars and the global-scale flows which are self-consistently established. The convection in these systems 
is richly time dependent and in our most rapidly rotating suns a striking pattern of spatially localized convection emerges. 
Convection near the equator in these systems is dominated by one or two patches of locally enhanced convection, with 
nearly quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong 
differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly 
rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation, as measured by 
latitudinal angular velocity contrast, Afl, increases with more rapid rotation while the relative shear, Afl/Q, decreases. 
In contrast, at more rapid rotation the meridional circulations decrease in both energy and peak velocities and break into 
multiple cells of circulation in both radius and latitude. 
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1 Rotation and Dynamo Action 

Rotation and magnetism couple with turbulent plasma mo- 
tions in stellar convection zones to drive dynamos and cy- 
cles of magnetic activity. When our sun was younger, it 
must have rotated more rapidly, as is suggested by obser- 
vations of rapidly rotating young solar-like stars and by the 
presence of the magnetized solar wind, which continually 
removes angular momentum from the star Observations of 
rapidly rotating stars indicate that generally more rapid ro- 
tation is correlated with stronger magnetism and perhaps 
stronger dynamo action. In the sun, global-scale dynamo 
action is thought to arise from the coupling of convection 
and rotation. Here we explore the effects of more rapid ro- 
tation on global-scale convective structures and in particular 
on the resulting differential rotation and meridional circula- 
tion. 

Helioseismology, which uses acoustic oscillations at the 
solar surface to probe the radial and latitudinal structure of 
the star as well as the convective flows beneath the sur- 
face, has revealed that the solar differential rotation pro- 
file observed at the surface prints throughout the convec- 
tion zone with two prominent regions of radial shear. The 
near-surface shear layer exists in the outer 5% of the sun, 
and the tachocline lies between the differentially-rotating 
convection zone and the deeper radiative interior, which is 



in nearly solid body rotation ( [Thompson et al ]l2003h . These 
shear layers, combined with the global-scale flows of merid- 
ional circulation and differential rotation, are thought to be 
the key components in the global solar dynamo which builds 
and rebuilds magnetic fields yielding 22-year cycl es of mag- 
netic activity. In the interf ace dynamo model (e . g. . ICharbonneau 



2005t IOssendrijveiil2003h . magnetic fields generated in the 



bulk of the convection zone are pumped downward into the 
tachocline where the strong radial shear builds large-scale 
fields that eventually erupt at the solar surface. The differ- 
ential rotation plays an important role both in building and 
organizing the global-scale fields while the meridional cir- 
culations may be important for returning flux to the base of 
the convection zone, enabling cycles of magnetic activity. 
The nature of internal differential rotation in other stars is 
unclear, though some observations of the surface differential 
rotation are now possible through a variety of technique s in- 



cluding photometric variability (Donahue et al. 1996; W alker et al 
Eoo7h . Dopple r imagin g (Donati et al. 2003) and Fourier trans- 
form methods (IReiners & Schmitt2003t) . We undertake sim- 
ulations to self-consistently establish differential rotation through 
the coupling of convection and rotation to study this crucial 
ingredient for stellar dynamo action. 
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2 Simulations of rapidly rotating suns 

To capture the essential couplings between convection, ro- 
tation and the resulting global-scale flows, we must employ 
a global model which simultaneously includes the spheri- 
cal shell geometry and admits the possibility of zonal jets, 
large eddy vortices and convective plumes which may span 
the depth of the convection zone. The solar convection zone 
is intensely turbulent and molecular values of viscosity and 
thermal diffusivity in the sun are estimated to be very small. 
As a consequence, numerical simulations cannot hope to re- 
solve all scales of motion present in real solar convection 
and a compromise must be struck between faithfully captur- 
ing the important dynamics within small regions and captur- 
ing the connectivity and geometry of the global scales. 

2.1 The ASH code 

Our tool for exploring stellar convection is the anelastic 
spher ical harmonic (AS H) cod e, which is described in de- 
tail in lClune et all (1 19991) and in lBrun et al.l h004 . ASH is a 
mature simulation code, designed to run on massively paral- 
lel architectures, which solves the three-dimensional MHD 
equations of motion under the anelastic approximation. The 
velocity field is fully nonlinear, but under the anelastic ap- 
proximation the thermodynamic variables are linearized about 
their spherically symmetric and evolving mean state with 
density p, pressure P, temperature T and specific entropy S 
all varying with radius. The overall density contrast is about 
25 between the top and bottom of the simulation. 

As a global-scale code, ASH cannot resolve all scales 
of motion present in real stellar convection zones. Instead, 
we conduct large eddy simulations (LES) which explicitly 
resolve the largest scales and describe scales below the grid 
resolution with subgrid-scale (SGS) modelling. Here we treat 
these scales with effective eddy diffusivities v and k, which 
represent momentum and heat transport by unresolved mo- 
tions in the simulations. For simplicity // and k are taken 
as functions of radius alone and in the simulations reported 
are proportional to p"^/^. The Prandtl number Pr = v/k 
is 0.25 for all simulations. 

2.2 Our simulations 

We are interested here primarily in the interaction between 
rotation and convection in deep stellar convection zones. In 
these simulations, we avoid the H and He ionization re- 
gions near the stellar surface as well as the stably strat- 
ified radiative interior and the tachocline of shear at the 
base of the convection zone. Thus we are concerned prin- 
cipally with the bulk of the convection zone and our com- 
putational domain extends from O.72i?0 to 0.96i?o. Solar 
values are taken for heat flux, mass and radius, and a per- 
fect gas is assumed. The reference or mean state of our ther- 
modynamic variable s is derived from a one-dimensional so- 
lar structure model (IBrun et al.ll2002h and is continuously 
updated with the spherically-symmetric components of the 



Table 1 Simulation parameter space. Evaluated at mid- 
depth are the Rayleigh number Ba, the Taylor number T^, 
the rms Reynolds number and fluctuating Reynolds 
number R'^, the Rossby number Ro and convective Rossby 
number Roc, and the diffusivities i/ and k in cm^/s. The ro- 
tation rate of the reference frame fto is given in multiples 
of Qq = 2.6 X 10~^ rad/s or 414 nHz. Also quoted are 
kinetic energies relative to the rotating system contained in 
convection (CKE), differential rotation (DRKE) and merid- 
ional circulations (MCKE) in units of 10® erg cm^"^ and as 
a percentage of the total. 



case 


Gl 


G5 




1^0 


5^0 


Ra 


3.22 X 10* 


1.29 X 10® 


Ta 


3.14 X 10^ 


6.70 X 10^ 


Re 


84 


543 


R'e 


63 


133 


Ro 


0.92 


0.28 


Roc 


0.61 


0.27 


V 


2.75 X 10^^ 


0.94 X 10^^ 


K 


11.0 X 10^2 


3.76 X 10^2 


CKE 


3.28 ( 59.0%) 


2.11(6.55%) 


DRKE 


2.26 ( 40.6%) 


30.1 ( 93.4%) 


MCKE 


0.025 ( 0.45%) 


0.007 ( 0.020%) 



thermodynamic fluctuations as the simulations proceed. We 
have conducted a number of 3-D hydrodynamic simulations 
for stars rotating from one to ten times the current solar rate. 
Here we concentrate on cases at the solar rate (case Gl) and 
at five times the solar rate (G5); the full suite of simulations 
are discussed at length inl Brown et al. ( 2007 ). The parame- 
ter space explored by these simulations is presented in Ta- 
ble [T] As we probe higher rotation rates, we decrease the 
eddy diffusivities in the simulation to maintain a compara- 
ble level of supercriticality. 



3 Convection in more rapidly rotating suns 

Examples of convective patterns in case Gl and G5 are illus- 
trated in Figure[T] The radial velocity near the top of the do- 
main is presented in Mollweide projection, which displays 
the entire layer with minimal distortion, with poles at top 
and bottom and the entire equatorial region displayed in the 
middle. The convection patterns are complex, time depen- 
dent and asymmetric. Owing to the density stratification, the 
narrow, fast downflow lanes are surrounded by broad, rela- 
tively weak upflows. Convection in the equatorial regions 
is dominated by large cells aligned with the rotation axis 
whereas at high latitudes the convection is more isotropic 
and cyclonic. With more rapid rotation, the convection cells 
decrease in horizontal size and align more strongly with the 
rotation axis but remain vigorous and turbulent. 

A striking finding in our more rapidly rotating systems 
is the pattern of spatially localized convection evident in 
Figure [T]:. This structure is a region of locally enhanced 
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Fig. 1 Evolution of global-scale convective flows with increasing rotation rate for case Gl (a,b) and case G5 (c,d). Shown 
in (a, c) are global views, in Mollweide projection, of radial velocity near the stellar surface (r — 0.95Rq). Upflows 
are light while downflows are dark, with scales indicated in m/s. At high rotation rates a striking pattern of modulated 
convection emerges at low latitudes, consisting of spatially modulated or patchy convection. Shown in {b, d) are azimuthal 
averages of angular velocity ft with radius and latitude. These have been further averaged in time over a period of roughly 
200 days. Plotted at right are radial cuts of angular velocity at selected latitudes as indicated. The black dashed contour in 
(d) denotes the constant propagation rate of the nests in such modulated convection. 



convection spanning from about ±30° in latitude. Here con- 
vective velocities and transport are enhanced, and strong 
plumes span the depth of the convection zone. This mod- 
ulated pattern emerges gradually with more rapid rotation 
and we have found states with either multiple patches of en- 
hanced convection or single isolated patches. At the high- 
est rotation rates, convection in the equatorial regions is 
present only within the patches, with complex streaming 
zonal flows in the interpatch regions. 

These active nests of convection persist for long periods 
of time (thousands of days) compared to the typical lifetime 
of individual convective cells (~ 10 — 40 days). They propa- 
gate in a prograde fashion relative to the stellar rotation rate 
and with a fixed angular velocity at all depths and latitudes 
in the equatorial region. The nests of convection are embed- 
ded in a strong radial shear, as shown in Figureflli. Near the 
surface the mean zonal flow streams through the nests while 
at depth the zonal flows are slower than the nests. The strong 
radial shear in the rapidly rotating cases leads to convective 
structures which are typically double-celled in radius. Near 
the stellar surface the individual convective cells propagate 
more rapidly than the nests, with these cells catching up and 
swimming through the structure. Near the base of the con- 



vection zone the nests propagate faster than the individual 
convective cells. 

4 Mean flows and circulations 

Differential rotation in our simulations is primarily estab- 
lished by Reynolds stresses from the convection. Convective 
enthalpy transport in latitude establishes a strong latitudi- 
nal temperature contrast in the more rapidly rotating stars, 
which are largely in thermal wind balance with their dif- 
ferential rotation. The profile of angular velocity realized in 
case Gl (Fig.[T})) is similar to the solar internal rotation pro- 
file, with a fast equator and a monotonic decrease in angular 
velocity towards the poles. The angular velocity profile re- 
alized in this simulation is significantly more columnar than 
the helioseismically deduced solar internal rotation profile, 
which is nearly cons tant on radial lines. Recent simulations 
dMiesch et al. 200 61) reveal that a more solar-like rotation 
profile can be achieved by including a latitude dependent 
thermal forcing that is consistent with a tachocline in ther- 
mal wind balance. As the nature of tachoclines and the mag- 
nitude of this thermal forcing is unclear in other stars, we 
have chosen here to neglect this latitude dependent thermal 
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Fig. 2 Evolution of circulations with faster rotation: (a) angular velocity contrast AQ /ileq plotted against bulk rotation 
rate H/Hq of the simulations in logarithmic scaling. A power law with exponent n = —0.6 is overplotted. {b) Peak 
meridional circulations at the top of the simulation (r = O.96i?0), showing a steep decline with more rapid rotation as the 
circulations break into multiple cells aUgned with the rotation axis. Here a power law with exponent n = —0.5 overlies the 
data. 



forcing and thus our bottom thermal boundary condition is 
one of constant flux. All of our rapidly rotating stars build 
angular velocity profiles with fast equators and slow poles. 
This strong differential rotation is realized even in the most 
rapidly rotating cases, which exhibit more extreme degrees 
of convective modulation. 



The kinetic energy contained in the differential rotation 
(DRKE) grows with rotation rate (Table [T]). This is accom- 
panied by a growth in the total latitudinal angular veloc- 
ity contrast A57 — ftcq — ^60° ■ This contrast grows more 
slowly than linearly with increasing rotation rate, and thus 
the relative contrast Aft/n decreases in the rapidly rotat- 
ing simulations (Fig.|2]2). Our simulations of rapidly rotat- 
ing suns appear to roughly follow a power law scaling at the 
highest rotation rates. However, the actual curve and scaling 
depends on the parameter space chosen for the simulations, 
and thus on the level of turbulence achieved, as discussed in 



detail in lBrown et alJ OOOTD 



5 Conclusions 

The complex coupling between rotation and convection drives 
a strong differential rotation in these more rapidly rotat- 
ing stars. The meridional circulations in contrast are much 
weaker, decreasing in both magnitude and relative contri- 
bution with more rapid rotation. As such, more rapidly ro- 
tating stars may have dynamos which differ substantially 
from the solar dynamo, where meridional circulations are 
thought to play an important role. These mean circulations 
are achieved despite the presence of strong spatial modu- 
lation in the convection of the equatorial regions. If these 
active nests of convection persist in stellar convection zones 
in the presence of magnetism, they may lead to features at 
the surface which propagate at rates distinct from that of 
the surface convection or the zonal flows of differential ro- 
tation. Such strong and persistent modulation in convective 
structures may have important observational consequences. 
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The energy contained in meridional circulations (MCKE) 
declines in both magnitude and percentage with more rapid 
rotation. The decrease in MCKE is accompanied by a de- 
crease in the magnitude of meridional velocities at the sur- 
face, as shown in Figure |2{), where we plot the magnitude 
of the peak time-averaged meridional velocity at the sur- 
face. This decrease in circulation strength is accompanied 
by a change in the topology of the circulations. In case Gl, 
the meridional circulations are largely single-celled in ra- 
dius and latitude in each hemisphere. In the more rapidly 
rotating stars, as in case G5, these cells break into multiple 
weak cells in both radius and latitude, and near the equator 
these cells align with the rotation axis. 
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